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TETHERED  BALLOON  MEASUREMENTS  MADE  AT 
NAVAL  AIR  DEVELOPMENT  CENTER 
(DECEMBER  1985) 


INTRODUCTION 

NADC  Scientists  have  built  a  remote  water  vapor  profiler. 
This  profiler  uses  a  laser  radar  (lidar)  to  measure  the  Raman 
shifted  signal  from  water  vapor  and  nitrogen.  The  ratio  of 
these  two  returns  provides  a  measure  of  the  mixing  ratio  allow¬ 
ing  this  parameter  to  be  profiled  from  near  the  surface  to  2 
km.  Since  the  lidar  operates  in  the  ultraviolet  (248  nm)  where 
ozone  absorption  is  a  maximum,  and  aerosol  scatter  is  signifi¬ 
cant,  the  ozone  and  aerosol  will  also  influence  these  lidar 
signals.  The  NADC  experiment  compared  the  lidar  water  vapor 
observations  with  standard  meterological  water  vapor  measure¬ 
ments. 


Front  5  through  11  December  1  935  ,  water  vapor  profiles 
were  observed  by  the  NADC  lidar,  NRL's  meterological  package, 
and  radiosondes  (operated  by  scientists  from  the  Pacific 
Missile  Test  Center,  Pt  Magu,  C  a ) .  This  report  discusses  the 
measurements  made  by  the  NRL  package  carried  aloft  by  a  53-ft. 
aerostat  balloon  from  the  surface  to  800  m.  (The  800  m  ceiling 
on  balloon  operations  result ed  from  F  A  A  restrictions  not 
balloon  limitations.)  Investigators  from  the  LTA  International 
Co.  operated  the  aerostat. 

The  instruments  attached  to  the  aerostat  measured  the 
water  vapor,  temperature,  ozone,  and  scattering  coefficient. 
Data  from  thermistors  (which  measured  the  wet  and  dry  bulb 
temperature  and  the  altimeter,  which  measured  pressure)  were 
combined  to  determine  the  water  vapor  mixing  ratio  .  In  addi¬ 
tion,  the  NRL  package  used  a  nephel oraeter  to  measure  the 
scattering  coefficient  and  a  sonde  to  measure  ozone  concentra¬ 
tions.  These  measurements  allowed  the  influence  of  aerosol 
scattering  and  ozone  on  the  lidar  operations  to  be  estimated. 
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INSTRUMENT  DESCRIPTION 


Reference  1  describes  the  NRL  instrument  package  in  de¬ 
tail.  Most  of  the  measurements  made  at  NADC  follow  the 
measurement  techniques  previously  discussed.  Only  the  ozone 
measurements  were  completely  new  for  the  NRL  package.  A 
modified  algorithm  with  a  temperature  correction  was  used  to 
determine  altitude  and  pressure.  Unlike  previous  experiments, 
the  freezing  of  the  'wet'  bulb  occurred  on  two  profiles  and,  on 
these  two  profiles,  the  water  vapor  mixing  ratio  was  determined 
from  the  ice  point  not  dew  point  relationships.  Only  a  brief 
instrument  description  is  included  in  this  report.  For  a  more 
detailed  description  of  the  instrument,  see  reference  1. 


ALTIMETER 

A  sensitive  barometric  altimeter  on  the  aerostat  balloon 
(Model  7000  Computer  Instrument  Corporation,  New  York,  N.  Y.) 
measures  instrument  height  as  the  difference  between  the 
pressure  at  the  instrument  altitude  and  the  surface  pressure. 
This  altimeter  follows  the  USSA-1962  standard  atmosphere  rela- 
ti onshi p  [2]  : 

p2  =  1013.3  (1-2.255X10-5  z2)5-256,  (1) 


where  p2  is  pressure  (mb)  at  height  z2  (m).  Near  the  surface, 
this  relationship  is  approximately: 


P2 


Pi  exp[-( f2-f 1)/k) 


(2) 


where  p 1  is  surface  pressure  (mb),  f2  is  altimeter  height  (ft), 
f^  is  surface  elevation,  and  k  is  a  constant  (2.464xlQ4).  The 
aerostat  elevation,  f2,  is  output  by  the  altimeter  and  directly 
converted  to  pressure.  The  constant  k  does  not  come  from  Eq  (1) 
but  from  a  calibration  done  by  Gerber;  however,  the  calibration 
results  are  consistant  with  Eq  (1).  The  temperature  and  water 
vapor  profiles  from  the  surface  to  the  altimeter  height  deter¬ 
mine  the  exact  pressure-height  relationship.  The  hypsometric 
equation  [3]  integrates  the  virtual  temperature  profile,  Tv,  to 
find  the  altimeter  height,  h: 


h(p2)-h(p1)  =  -f  2  R  Tv(p)  g"1  d  ln(p), 
;Pl 


(3 


where  R  is  the  dry  gas  constant,  and  g  is  the  gravitional 
constant.  Cool  conditions  during  the  NADC  experiment  made  the 
water  vapor's  contribution  to  virtual  potential  temperature 
negligible;  thus,  the  temperature  profile  was  used  instead  o- 
the  virtual  temperature  profile  changing  the  pressure-height 
relationship  to: 


R  T  (  f )  g"  i  k"1  df  , 


(4) 


h(  f  2 ) 


h(  f  ! ) 
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where  T  is  the  temperature  (K).  Eqs.  (2)  and  (4)  were  tested 
on  flights  8a  (12/8/85)  and  12b  (12/10/85)  when  the  boundary 
layer  appeared  to  be  well  mixed.  The  lapse  rate  in  these 
temperature  profiles  was  extremely  close  to  adiabatic.  On 
future  experiments,  the  altimeter  should  be  recalibrated  to 
ensure  the  'constant'  has  not  drifted. 


PSYCHROMETER 

The  wet  bulb  temperature,  dry  bulb  temperature,  and  pres¬ 
sure  measurements  determine  the  water  vapor  mixing  ratio.  The 
altimeter  readings,  Eq.  (4),  are  converted  to  pressure.  Matched 
thermistors  mounted  inside  two  concentric  heat  shielding  tubes 
measure  the  wet  and  dry  bulb  temperatures.  By  being  in  line 
with  the  aerostat  balloon,  these  tubes  allow  the  wet  bulb  to  be 
wind  aspirated.  The  water  vapor  mixing  ratio,  W  (gm  water 
vapor/kg  dry  air),  is  derived  from: 

W  =  103  (Mv/Md)  (e/(p-e))  ,  (5) 

where  (Mv/Md)  is  the  ratio  of  molecular  weights  of  water  and 
dry  air  (0.622),  e  is  the  atmospheric  partial  pressure  of  water 
vapor.  The  relationship  [4]  estimates  the  water  vapor 
pressure,  eg  (mb),  at  temperature  T  (K)  by  using  : 

1 og(e0)  =  -7.90298(Ts/T-l)+5. 02808  log(Ts/T) 

-  1.3816X10“7[1011-334  (I-T/TsJ.j] 

+8.1328xl0-3[10“3*49149  (Ts/T-1)-1]  +  log(eos)  ,  (6) 

where  T  s  is  the  steam  point  (373.16  K),  eos  is  saturation  vapor 
pressure  at  the  steam  point  (1013.246  mb),  and  log  is  the  loga¬ 
rithm  base  10.  When  the  'wet'  bulb  is  frozen,  the  vapor 
pressure  is  calculated  from  the  ice  point: 

1 og ( eg )  =  -9.09718  (Tg/T-1  )  -3.56654  1  o g ( T 0 / T )  + 

0.876793  (1-  T/Tg)  +  log(ei0),  (7) 

where  e^0  is  the  saturation  at  ice  point  (6.1071  mb)  and  Tg  is 
the  ice  point  (273.16  K).  The  vapor  pressure,  eg,  observed  by 
the  wet  bulb  and  determined  by  either  Eq  (6)  or  (7),  is 
correct ed  for  ventilation  effects: 

e  «  e«  -  [0.00066  (  1  +  0.001  1  5  T,.,)l  ( TH  -  T ... )  p  ,  (8) 


where  and  Tw  are  the  temperatures  (K)  read  by  the  dry  bulb 
and  wet  bulb  respectively.  When  the  vapor  pressure  is  calcu¬ 
lated  from  a  •frozen  bulb,  the  constants  in  eqn  (8)  are 
multiplied  by  C-882,  the  ratio  of  latent  heat  of  evaporation  of 
water  to  that  o'  ice. 

0Z0NES0NDE 

A  lightweight  electrochemical  concentration  cell  in  a 
sonde  hung  below  the  aerostat  creates  a  current,  which  is 
proportional  to  the  ozone  partial  pressure,  P3  (mb): 

p3  =  4 .307  xlO"3  I  Tb  t,  (9) 


where  I  is  the  sensor  current  (micro-amps),  Tb  is  sonde 
temperature  (K),  and  t  is  time  (s)  the  pump  takes  to  force  100 
ml  of  air  through  the  sonde  (approximately  29.5).  For  plotting, 
the  ozone  partial  pressure  is  divided  by  the  local  pressure  and 
converted  to  parts  per  billion  (ppb).  Prior  to  each  flight,  the 
sonde  received  new  chemical  solutions.  A  known  ozone  concentra¬ 
tion  is  pumped  through  the  sonde  and  the  sonde  reading  for  this 
concentration  is  tested.  After  the  ozonesonde  is  mounted  on  the 
aerostat  and  before  the  balloon  is  launched,  the  sonde  is  run 
for  five  minutes  with  and  five  minutes  without  an  ozone  free 
filter  to  determine  the  sonde's  zero  reading. 

NEPHELOMETER 

The  nephel ometer  uses  a  photographic  flash  lamp  to 
uniformly  illuminate  the  inside  of  an  integrating  sphere.  The 
scattered  light  is  measured  by  a  photomultiplier,  which  has  a 
field  of  view  that  goes  through  the  center  of  the  integrating 
sphere  and  is  carefully  aligned  with  a  light  trap  on  the  oppo¬ 
site  side  of  the  sphere,  ensuring  that  no  directly  reflected 
light  enters  the  photomultiplier.  A  photopic  interference 
filter  on  the  photomultiplier  allows  the  nephelometer  to 
measure  visibility  (using  Koshmieder's  critera).  The  scattering 
coefficient,  b$  (l/km),  measured  by  the  nephelometer  is: 

b$  =  1.22  (Vout  -  0.107),  (10) 

where  V  Q  u  ^  is  the  output  voltage.  While  the  nephelometer  can 
make  readings  when  visibility  ranges  from  0.330  to  90  km,  the 
scattering  volume  is  small  (  1  cm-)  and,  when  visibility  is 
greater  than  about  50  km,  the  statistical  properties  of  the 
scattering  can  vary  in  short  time  periods.  To  minimize  the 
variance,  the  measurements  were  avera ged  for  2  minutes. 


ANEMOMETER 


A  cup  anemometer,  mounted  on  a  tail  fin  of  the  aerostat 
balloon,  transmits  pulse  counts,  which  have  been  calibrated  to 
yield  the  wind  speed,  u: 

0.357  P  -  7.86  P > 36 
u  =  ,  0 

3 .86X10"-5  P  <  36  ,  (11 

where  P  is  the  pulse  count.  The  wind  direction  could  be  esti¬ 
mated  from  the  orientation  of  the  aerostat  balloon;  however, 
wind  direction  was  not  an  objective  of  the  experiment  and  was, 
therefore,  not  recorded.  Flights  6  and  8  showed  more  wind  speed 
variation  than  expected  for  the  meterol ogi cal  conditions.  This 
variance  could  be  increased  by  the  aerostat  motion  and  balloon 
sheltering  of  the  anemometer. 


OBSERVATIONS 


Aerostat  flights  were  conducted  on  every  day  from  5  to  11 
December  1985  except  on  6  December  when  high  winds  kept  the 
aerostat  from  flying.  There  were  two  flights  on  all  flight 
days,  except  the  last,  one  in  the  morning  and  one  in  the  after¬ 
noon.  Flight  1  was  used  to  test  instrumentation.  Altimeter 
problems  restricted  measurements  on  Flight  2.  The  next  eleven 
flights  which  produced  useful  profiles  are  described  in  Table  1 
and  will  be  discussed  in  more  detail  later  in  this  report  (for 
figures  see  Appendix  A).  After  the  first  two  flights,  the  only 
restrictions  came  from  high  winds  and  airport  operations.  When 
time  permitted,  profiles  were  made  to  800  m. 

The  aerostat  was  launched  from  the  apron  of  the  NADC 
runway  (40°)  12'  30"  N  and  75°  5'  E,  see  Fig.  1)  which  was  114 
m  above  sea  level.  The  two-  to  three-story  NADC  headquarters 
building  was  110  m  northwest  of  the  launch  site  and  a  hangar 
facility  was  80  m  northeast  of  the  launch  site.  The  runway  was 
southeast  of  the  launch  site  with  the  buildings  800  m  from  the 
site.  Pennsylvania  State  Highway  132  was  southwest  (350  m) 
from  the  experimental  site  and  the  nearest  buildings  were  580  m 
from  the  experimental  site.  Outside  NADC,  the  land  has  been 
developed  as  residential  and  light  industry.  The  Pennsylvania 
State  Turn  Pike  (1-70)  was  4  km  to  the  south. 

The  aerostat  made  c o n s t a n t - r a t e  ascent  and  descent  pro¬ 
files  and  was  held  at  constant  altitudes  to  determine  the 
water  vapor  mixing  ratio  statistics.  Those  statistics  are 
listed  in  Table  2.  Profiles  of  the  water  vapor  mixing  ratio, 
ozone,  scattering  coefficient,  wind  speed,  virtual  potential 
temperature,  and  natural  logarithm  of  the  scattering  coeffi¬ 
cient  are  plotted.  The  plots  are  labeled  with  the  flight  number 
and  letters  denoting  the  separate  ascent/d escent  profiles.  Most 
fights  involved  ascents  to  a  maximum  altitude  followed  by 
descent  to  the  surface.  On  these  flights,  the  ascent  and 
descent  were  marked  'a'  and  *  b '  respectively.  Flights  9,  11, 
and  13  involved  direct  ascents  to  maximum  altitude,  a  descent 
to  intermediate  altitudes,  ‘re-ascent*  to  the  maximum  altitude, 
and  a  final  descent  to  the  surface;  these  profiles  are  labeled 
‘a1,  ‘ b '  ,  '  c '  ,  and  ' d  '  respectively.  The  aerostat  heights  as  a 
function  of  time  are  plotted  allowing  easy  visualization  of 
balloon  operations.  The  wind  speed  profiles  were  not  a  major 
objective  of  this  experiment  and  are  included  for  informational 
reasons . 


In  general,  the  flights  were  uneventful  and  the  profiles 
of  water  vapor  mixing  ratio  have  very  little  structure.  Flight 
13  was  a  major  exception  with  an  increase  of  the  mixing  ratio 
from  4.5  to  6.5  gm/kg  at  300  m.  A  rapid  decease  from  3  to  2 
gm/kg  was  observed  at  the  top  of  Flight  7.  The  top  of  the  3 
gm/kg  layer  deceased  during  the  profile. 
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Table  1  Flight 

Summary 

Time 

Fixed  Heights 

(EST  ) 

(m  ) 

0723-0945 

334,454,570,691,786 

693,573,457,333 

1220-1432 

344,463,583,702,810, 

708,590,475 

0745-0942 

287,391,487,585,691 

1217-1436 

309,412,533,640,726 

0753-0928 

730,625 

1251-1431 

761,617,540 

0731-0943 

525,760 

1430-1530 

741 

0747-0943 

746,321,523,743 

1203-1258 

637 

0735-0945 

759,304,598,771 

EXPERIMENTAL  SITE 


TABLE  2  -  FIXED  LEVEL  STATISTICS 


Time 

(HHMM  ) 

Height 

(m  ) 

Mixing  Ratio 

(gm/K) 

tart 

Stop 

Mean 

Va  r 

Mean 

Va  r 

12/05/85 

741 

0749 

334.0 

0.6 

1.43 

0.01 

753 

0807 

454.1 

0.5 

1.27 

0.03 

808 

0818 

570.3 

0.8 

1.11 

0.13 

822 

0831 

691.5 

0.9 

0.90 

0.06 

835 

0845 

7  86.1 

0.8 

0.76 

0.05 

848 

0858 

693.6 

1.2 

1.11 

0 .14 

1902 

0912 

573.3 

1.4 

1.27 

0.04* 

i916 

0926 

456.8 

0.9 

1.70 

0.04* 

1930 

0940 

333.1 

0.5 

2.09 

0.06* 

.233 

1243 

343.8 

0.5 

1.47 

0.10 

.247 

1256 

463.1 

0.4 

1.25 

0.02 

.301 

1310 

583.0 

0.8 

1.22 

0.19 

.314 

1322 

701.7 

0.6 

0.94 

0.01 

.328 

1336 

810.4 

1.2 

0.59 

0.01 

'339 

1358 

707.5 

0.7 

0.86 

0.07 

L402 

1412 

590.1 

0.6 

1.23 

0.03 

[416 

1428 

475.0 

0.6 

1.18 

0.03 

12/0  7/85 

3757 

0806 

287.4 

1.6 

2.39 

0.01 

3811 

0821 

391.3 

2.6 

2.34 

0.01 

3826 

0835 

487.0 

3.2 

2.25 

0.01 

0840 

0848 

584.7 

4.2 

2.24 

0.03 

0854 

0912 

690.6 

7.4 

2.11 

0.06 

1232 

1247 

308.5 

1.9 

2.46 

0.02 

1250 

1305 

412.1 

4.8 

2.43 

0.01 

1310 

1318 

53  2.9 

6.0 

2.44 

0.01 

1327 

1342 

639.7 

3.4 

2.48 

0.01 

1350 

1404 

725.7 

4.1 

2.50 

0.01 

♦Frozen  'Wet'  Bulb 


TABLE  2  (Cont’d)  -  FIXED  LEVEL  STATISTICS 

Time 

(HHMM  ) 

Height  (m  ) 

Mixing  Ratio 

(gm  /K) 

tart 

Stop 

Mean  Var 

Mean 

Va  r 

12/08/85 

816 

0838 

730.3  14.0 

2.31 

0.09 

850 

0904 

625.4  4.6 

2.43 

0.23 

320 

1341 

760.5  10.3 

2.74 

0.03 

345 

1356 

616.9  6.8 

2.92 

0.03 

410 

1416 

539.8  6.2 

2.94 

0.03 

12/09/85 

813 

0847 

524.6  0.8 

3.22 

0.03 

859 

0914 

759.5  1.2 

3.52 

0.01 

456 

1514 

740.7  34.5 

3.37 

0.11 

12/10/85 

316 

0824 

746.4  1.6 

3.37 

0.02 

838 

0847 

320.5  0.5 

3.41 

0.01 

354 

0905 

523.2  0.6 

3.50 

0.02 

913 

0917 

742.9  0.4 

3.30 

0.01 

229 

1236 

636.8  8.0 

3.46 

0.05 

12/11/35 

hKM 

Itflil 
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Measurement  problems  were  restricted  to  Flights  3,  5,  and 
11.  The  wet  bulb  froze  during  the  descents  of  Flights  3  and  5 
at  heights  573  and  690  m  respectively.  The  gradual  increase  in 
mixing  ratio  during  the  descent  on  these  flights  may  be  caused 
by  drying  of  the  'ice'  bulb  by  sublimation.  On  flight  11,  the 
wet-bulb  water  container  was  dry  until  just  before  launch,  and 
the  wet-bulb  equilibrium  at  the  maximum  altitude;  therefore, 
the  water  vapor  profile  11a  is  omitted. 

The  ozonesondes  worked  well  throughout  the  observation 
period  and  seemed  tied  to  boundary  layer  dynamics.  For  example, 
on  Flight  5,  the  ascent  showed  a  layer  in  which  the  ozone  had 
apparently  been  depleted  overnight,  but  by  the  descent  ozone 
was  added  to  this  layer  by  mixing.  While  the  ozone  concentra¬ 
tions  fluctuated  wildly  on  Flights  3  and  4,  the  measurements 
appear  accurate.  The  two  flights  were  made  with  different 
ozonesondes.  (NOAA  loaned  us  two  sondes).  The  calibrations 
prior  to  these  flights  were  normal  and  the  large  variations 
appeared  only  in  the  ozone  measurements.  The  ozonesonde 
temperature  changed  slowly  with  altitude.  The  only  problem  in 
the  ozone  measurements  was  anomalous  negative  ozone  readings  at 
the  surface.  On  Flights  7,  9,  11,  12,  and  13,  the  surface  ozone 
readings  with  the  ozone  free  filter  were  higher  than  the  sur¬ 
face  reading  without  the  filter.  These  negative  readings  were 
probably  the  result  of  sulfur  dioxide  con  t  am  i  na  t  i  on  [5].  The 
negative  ozone  readings  were  extremely  small  and  are  not 
expected  to  significantly  influence  the  ozone  profiles. 

The  time-height  relationship  of  the  aerostat  should  be 
taken  into  account  when  the  profiles  are  analyzed.  The  ascent 
rate  and  frequent  stops  to  measure  statistics  at  constant 
heights  change  the  appearance  of  the  profiles.  On  Flight  3,  the 
ascent  to  maximum  altitude  took  1.5  hours.  An  'inversion' 
appears  on  profile  6a;  however,  this  inversion  was  caused  in 
part  by  the  warming  that  took  place  while  the  aerostat  was  held 
at  300  m . 
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SUMMARY 


The  NRl  package  performed  almost  flawlessly  making  mea¬ 
surements  for  eleven  flights  and  over  22  flight  hours  without  a 
major  instrumental  failure.  On  many  days,  the  water-vapor 
mixing  ratio  was  almost  constant.  Data  comparisons  on  these 
days  should  be  possible  without  the  spatial  and  temporal  sepa¬ 
rations  between  instruments  introducing  large  spurious 
differences  between  the  measurements.  Profiles  from  Flights  7 
and  13  have  'step  functions’  that  should  be  apparent  on  all 
three  instruments. 


The  ozonesonde  required  considerably  more  work  than  the 
other  instruments  on  the  NRL  package;  however,  the  results 
amply  justified  this  effort.  Ozone  was  frequently  the  best 
indicator  of  the  boundary  layer  and  boundary  layer  mixing.  The 
apparent  discovery  of  sulfur  dioxide  (if  confirmed  by  other 
observations)  may  have  significance  beyond  contaminating  ozone 
measurements.  As  a  tri -atomic  molecule,  it  absorbs  in  the 
ultra-violet  and  could  also  influence  lidar  operations.  The 
purchase  of  an  ozone  detector  as  an  addition  to  the 
meteorological  package  should  be  seriously  consider ed. 

This  experiment  illustrates  the  importance  of  the  NRL 
meteorological  package  and  the  aerostat  balloon.  When  used  with 
a  real-time  display,  the  package  hung  below  a  tethered  balloon 
has  the  ability  to  map  m e t e r ol o g i c a  1  and  optical  parameters. 
Accurate  measurement  of  these  parameters  is  essential  for  both 
future  boundary  layer  experiments  and  optical  instruments 
tests . 
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APPENDIX  A  -  Profiles  of  Mixing  Ratio,  Ozone,  Scattering  Ratio 
Virtual  Potential  Temperature,  and  Wind  Speed 

The  measurements  made  by  the  'Aerostat'  tethered  balloon 
are  plotted  as  a  function  of  height  for  the  flights.  The  pro¬ 
files  are  grouped  by  flights  in  chronological  order  from  Flight 
3  to  13.  Each  plot  has  a  flight  number,  date,  and  time  in  the 
upper  right  corner.  The  first  plot  for  each  flight  is  a  time- 
height  profile  followed  by  the  mixing  ratio,  ozone,  scattering 
ratio,  virtual  potential  temperature,  and  wind  speed  profiles. 
The  scattering  ratio  has  a  natural  logarithm  scale.  All  other 
plots  have  linear  scales. 
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